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ABSTRACT 

Swift observations suggest that the central compact objects of some gamma-ray bursts (GRBs) could be newly born 
millisecond magnetars. Therefore, by considering the spin evolution of the magnetars against r-mode instability, we 
investigate the role of the magnetars in GRB X-ray afterglow emission. Besides modifying the conventional energy 
injection model, we pay particular attention to the internal X-ray afterglow emission, whose luminosity is assumed 
to track the magnetic dipole luminosity of the magentars with a certain fraction. Following a comparison between 
the model and some selected observational samples, we suggest that some so-called "canonical" X-ray afterglows 
including the shallow decay, normal decay, and steeper-than-normal decay phases could be internally produced by the 
magnetars (possibly through some internal dissipations of the magnetar winds) , while the (energized) external shocks 
are associated with another type of X-ray afterglows. If this is true, from those internal X-ray afterglows, we can 
further determine the magnetic field strengths and the initial spin periods of the corresponding magnetars. 
Subject headings: gamma ray burst: general — stars: neutron 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are short, intense flashes of 
soft gamma-rays (~ 0.01 — 1 MeV), which are always fol- 
lowed by long-lasting low-frequency afterglow emission. 
Usually, the afterglow emission is attributed to an ex- 
ternal forward shock arising from the interaction of the 
GRB outflow with the circumburst medium, whereas the 
mechanisms responsible for the bursts are under more de- 
bate. Since the launch of the Swift spacecraft (Gehrels 
2004), many detailed features of the GRB X-ray after- 
glows have been revealed by the X-Ray Telescope (XRT) 
aboard. Then Nousek et al. (2006) and Zhang et al. 
(2006) phenomenologically summarized a "canonical" X- 
ray light curve (LC) with four smooth segments (some- 
times superposed by some sharp flares), although the 
LCs owning all components are actually in the tiny mi- 
nority. Specifically, the four different emission phases de- 
fined by them include: (1) Initial steep decay phase that 
is widely accepted to be the tail of the prompt emission 
(i.e., curvature effect; Fenimore et al. 1996; Kumar & 
Panaitescu 2000). (2) Shallow decay (even a plateau) 
phase that is usually ascribed to a continuous energy in- 
jection into the external shock (Rees & Meszaros 1998; 
Dai k Lu 1998a,b; Zhang & Meszaros 2001). (3) Normal 
decay phase that does not contradict with the standard 
external shock model. (4) Steeper-than-normal decay 
phase that is often connected to the jet break (Rhoads 
1997). To summarize, in the conventional picture as de- 
scribed above, the emission during the shallow decay, 
normal decay, and steeper-than-normal decay phases (of 
interest in this paper) are usually considered to be as- 
sociated with the external shock, while the initial steep 
decay emission as well as the flares are probably of inter- 
nal origin. 

The X-ray shallow-decay and flare emission strongly 
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suggest that GRB central objects should have long ac- 
tivities after the bursts. Therefore, highly magnetized, 
rapidly spinning pulsars (i.e., millisecond magnetars) 
gradually become a popular candidate of the central com- 
pact objects (e.g., De Pasquale et al. 2007; Metzger et 
al. 2007; Troja et al. 2007; Zhang & Dai 2008, 2009; 
Bucciantini et al. 2009; Corsi & Meszaros 2009; Lyons 
et al. 2009), although the black hole model is still an 
attractive choice. It is a difficult task to observationally 
distinguish between the magnetar and black hole mod- 
els. However, some magnetohydrodynamic simulations 
showed that the magnetar model could be more mature 
in the sense that it provides quantitative explanations for 
the durations, energies, Lorentz factors, and collimation 
of long GRB outflows (Metzger 2010). Moreover, with a 
spinning-down magnetar, some Swift-XKF features can 
be explained well (e.g., Dai et al. 2006; Fan & Xu 2006; 
Yu & Dai 2007) . Especially, the unusual X-ray after- 
glow LC of GRB 070110, where a nearly constant X-ray 
emission is followed by a very steep decline of a ~ 9 (a 
is the decline index of t~ a ), can be understood by as- 
cribing the plateau emission to magnetar- driven internal 
emission (Troja et al. 2007). 

The internal plateau emission of GRB 070110 also tells 
us that, in addition to the initial steep decay and flare 
phases, the X-ray emission during any other afterglow 
phases could also arise from some internal dissipation 
mechanisms. Therefore, it is fair to consider that the 
long active central objects of some other GRBs could also 
play an essential, relatively more direct role in their af- 
terglow emission, at least in the X-ray band. In some ex- 
treme situations, we suspect that such an internal-origin 
emission component could even dominate the total X-ray 
afterglow emission of a GRB, whereas the external shock 
emission is outshined in the X-ray band. In other words, 
the conventional external shock model is only one choice 
among various afterglow origin models, as also proposed 
by some authors before (Ghisellini et al. 2007; Kumar 
et al. 2008; Cannizzo & Gehrels 2009; Lindner et al. 
2009; Lyutikov 2009). A wide investigation on the GRB 
X-ray afterglows given by Willingale el al. (2007) indeed 
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showed that more than one hundred X-ray afterglows can 
be divided into two emission components, one of which 
is probably of internal origin. 

Following the observational results and the above two 
theoretical considerations, in this paper, we investigate 
in more detail the role of millisecond magnetars in the 
X-ray afterglow emission of some GRBs, based on a care- 
ful analysis on the spin evolution of the magnetars. Be- 
sides modifying the conventional energy injection model 
(Dai & Lu 1998a, b; Zhang & Meszaros 2001), we pay 
particular attention to the internal X-ray afterglow emis- 
sion that is produced by the magnetars, possibly through 
some internal dissipations in the magnetar winds. The 
luminosity of this internal emission is assumed to simply 
track the magnetic dipole luminosity of the magnetars 
with a certain fraction. To summarize, the observed X- 
ray afterglows could be emitted from two different regions 
(i.e., an internally-dissipated magnetar wind and an en- 
ergized external shock) at very different radii. The com- 
petition between these two emission components leads to 
a diversity of the X-ray afterglow LCs. 

In Section 2, we briefly review the spin-down of magne- 
tars against r-modc instability. We analyze the temporal 
behaviors of X-ray afterglows by combining the contri- 
butions from a magnetar wind and an external shock in 
Section 3, where the energy injection from the wind to 
the shock is also taken into account. In Section 4, some 
observational samples arc selected and fitted in order to 
confront the model with observations. Meanwhile, some 
implications to the magnetars from the afterglow data 
are discussed. Finally, a summary and discussion are 
given in Section 5. 

2. SPIN-DOWN OF MAGNETARS 

For long GRBs associated with Type Ibc super- 
novae, the central magnetars could be formed during the 
core-collapse of massive, rotating stars, since powerful 
magneto-centrifugal outflows from the nascent magne- 
tars may stave off black hole formation entirely (Met- 
zger 2010). For short GRBs, the merger of compact bi- 
naries may also give rise to a massive neutron star if 
the equation of state of stellar matter is stiff enough, as 
implied by the observations of kilohertz quasi-periodic 
oscillations in accreting neutron stars (Kluzniak 1997; 
Kluzniak & Ruderman 1998). So it seems somewhat ac- 
ceptable that the central compact objects of some GRBs 
are millisecond magnetars. The idea that GRBs could 
originate from magnetars actually had been proposed 
almost two decades ago (Usov 1992, 1994; Paczynski 
1992; Duncan & Thompson 1992), far before the Swift 
era. Specifically, several central engine mechanisms have 
been proposed, e.g., a neutrino-driven wind (Thompson 
1994; Metzger et al. 2007; Bucciantini et al. 2009;), a 
magnetic reconnection-accelerated wind (Drcnkhahn & 
Spruit 2002), or an hyperaccretion onto the magnetars 
(Zhang & Dai 2008, 2009). Anyway, these violent pro- 
cesses lead the initial spin evolution of the magnetars 
during the first tens of seconds to be very complicated. 

However, on the relatively longer (afterglow) 
timescales concerned here, the short-term processes 
as mentioned above could be no longer important. Then 
the spin-down of the magnetars would be mainly due to 
electromagnetic torque and the torque connected with 
gravitational wave radiation. For the latter, first, a 



magnetic-field-caused equatorial cllipticity can generate 
gravitational-quadrupole radiation, which however can 
not provide an essential contribution to the spin-down 
unless with an extremely high magnetic field (Usov 
1992). Secondly, a relatively stronger gravitational wave 
radiation can be produced through some nonaxisymmet- 
ric stellar perturbations. For nascent neutron stars, such 
perturbations can be easily created by r-mode instability, 
which arises from the action of the Coriolis force with 
positive feedback (Andersson 1998; Friedman & Morsink 
1998) succumbing to gravitational radiation-driven 
Chandrasekhar-Friedman-Schutz instability. 

Following a phcnomenological second-order model for 
the r-mode evolution (Owen et al. 1998; Sa 2004), we 
can calculate the spin evolution of a magnetar by (Sa 
2004; Yu et al. 2009a) 

4a 2 P P 
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where P is the spin period of the magnetar, a is the 
dimcnsionless amplitude of the r- modes, and 5 is a free 
parameter describing the initial degree of the differen- 
tial rotation of the star. The gravitational and mag- 
netic braking timescales can be written as r g = 144P^ 3 

s and T m — 4 x 10 5 B^ 2 P 2 3 s, respectively 3 . The viscous 
damping of the r-modes is ignored here due to the high 
temperatures of ~ 10 10 K in nascent neutron stars. Cor- 
respondingly, the evolution of the r-mode amplitude can 
be calculated from (Sa 2004; Yu ct al. 2009a) 

2a 2 1 
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In the case of T g j » r m .i (i.e., B > B c = 5 x 10 5 P 4 _ 3 
G), the spin-down would be dominated exclusively by the 
magnetic dipole radiation and the r-modes cannot arise 
sufficiently rapidly. In this case, the spin evolution can 
be easily expressed as usual as P(t) = Pj(l + t/Tm) 1 / 2 
with T m = 2x 10 5 B^P?_ 3 s. 

In contrast, for not very high magnetic fields (i.e., 
Tg.i <C T m ,i), the spin-down should be first dominated 
by the gravitational wave radiation. By ignoring the 
magnetic term, Equations. (1) and (2) can be solved 
analytically and an asymptotic solution can be written 
as (Sa & Tome 2005, 2006) 



P(t) 



Pi [1- ^a 2 l {S + 2)cxp{2t/T g , l )] \ fort<T 3 , 
1.6P(Vr 9 . 2 ) 1/5 , for t > T g . 

, , (3) 

Here the break time can be solved from (d a/dt ) t= T g — 
to be T g = ax 10 3 P l 6 _ 3 s, where the prefactor 0.7 < a < 
3.2 for a wide parameter region of 10~ 10 < a < 10~ 6 and 
< S < 10 8 (Yu et al. 2009b). As the increasing of the 
spin period, the magnetic braking effect would eventually 
exceed the gravitational braking effect and then the time- 
dependence of the spin period changes from P oc i 1 / 5 to 
P oc i 1//2 . We denote the change time by T c , which will 
be given explicitly later. 

3 Hereafter, some basic structural parameters such as the mass, 
the radius, and the moment of inertial of the magnetar are taken 
to be 1.4Mq, 10 6 cm, and 10 45 g cm 2 , respectively. Additionally, 
the convention Q x = Q/10 x is adopted in cgs units and a subscript 
"i" represents the initial values of the quantities. 
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3. GRB X-RAY AFTERGLOW EMISSION 

After the bursts, the central remanent magnetars could 
still keep abundant rotational energy, although a consid- 
erable fraction has been expanded on the bursts. This 
remaining energy could be released peacefully and per- 
sistently, and drive a continuous magnetar winds. We 
believe that the magnetar winds are probably able to 
produce long-lasting emission (i.e., afterglow emission), 
through some internal dissipation mechanisms such as 
magnetic reconnection (e.g., Giannios & Spruit 2005) 
and the terminative shock of the winds (e.g., Dai 2004; 
Yu & Dai 2007) etc. Such an internal-origin emission 
component would of course compete with the external 
shock emission component, while the latter alone is usu- 
ally challenged by some multi-wavelength afterglow ob- 
servations (e.g., chromatic LC breaks). 

In order to compare with X-ray observations, we first 
simply assume that the (isotropically-equivalent) lumi- 
nosity of the wind X-ray afterglow emission tracks the 
magnetic dipole luminosity of the magnetars with a con- 
stant fraction (£; X-ray radiation efficiency) as 

£T = £W/u = 10 47 e-i/ B ,-i L md,47 erg s-\ (4) 

where fs = (1 — cos6 w ) is the beaming factor of the 
magnetar winds with 9 W being the half-opening angle of 
the winds. To be specific, the magnetic dipole luminosity 
of a magnetar can be calculated by 

TV 2 

imd(t) = 10 47 F(t)B 2 14 Pr* 3 erg a" 1 , (5) 

where / is the moment of inertial of the star and Q = 
2n/P is the spin frequency. According to the analysis 
on the spin-down of the magnetars, the time-dependence 
F(t) of the magnetic dipole luminosity can be expressed 
approximately as 
For B > B c , 




For B < B c , 



F(t) 



t <T„ 
t>T„ 



t<T g , 

Tg<t< T c , 

t>T c , 



(6) 



(7) 



where T c = (T^/T^j) 1 ^ 2 '^ . Following Equation ©, the 
value of q can be taken as 0.8 approximately (a more ex- 
act numerical calculation would show q ~ 1). As shown 
by Equation (J7J for B < B c , we surprisingly find that 
the double-broken power-law behavior of C 1 ^" is in good 
agreement with the so-called canonical X-ray LC includ- 
ing all of the shallow decay, normal decay, and steeper- 
than-normal decay phases. This may be just a coinci- 
dence, but alternatively it seems also acceptable to as- 
cribe some observed double-broken power-law X-ray af- 
terglow emission to the magnetar winds rather than the 
external shocks. In this case, the change in the slopes of 
the LCs is just because of the evolution of the magnetic 
dipole luminosity, but independent of the external shock 
physics. As an analogical consideration for the magne- 
tars with B > B c , Equation ([6]) also predicts some X-ray 



LCs behaving as a plateau followed by a steep decay with 
a ~ 2. However, the rare observation of such a type of 
X-ray LCs implies that the magnetic field strengths of 
the GRB magnetars are generally lower than the criti- 
cal value B c . 4 So we would mainly concern the case of 
B < B c in the following calculations. 

Although the internal X-ray afterglow emission is sug- 
gested above, we still think that the external shock could 
play an important role in the afterglow emission. In the 
calculation of shock dynamics, we also take into account 
the energy injection into the shock from the wind as 
considered before (e.g., Dai & Lu 1998a, b; Zhang & 
Meszaros 2001), even though a part of the wind energy 
could have been radiated directly. Following Equations 
([5]) and (J7J for B < B c , the increasing (isotropically- 
equivalent) shock energy can be approximately written 
as (for q < 1) 



£{t) »£jX< 



i m ft 



1-9 
1-9 



t < T el , 

T C1 <t< Tg, 



T g <t<T c 



(8) 



£f, t>T c 



The starting time of the energy increase can be calculated 
by T ci = Si/iCfs 1 ^) = lOOCo^.so/s.-i^dM? s ' 
where £ < (1 — £) is the fraction of the wind energy that 
is injected into the shock. Equation © shows that the 
energy injection process can be divided into two stages 
for q < 1. However, since the value of q is actually very 
close to one, the energy increase of the shock after T g 
is difficult to be detected and thus can be neglected. 
Therefore, in comparison with the conventional energy 
injection model proposed by Dai & Lu (1998a, b) and 
Zhang & Meszaros (2001), the main difference here is 
that the duration of the energy injection is determined 
by T g rather than T m . For shock dynamic evolution 5 , 
we first introduce the deceleration time of the shock as 
Tdoc = ^l5o^i,2.3 n o 1 s (< T a < T g)- Before T doc , 
the shock deceleration can be neglected (i.e., T ~ Ti) 
since A4 SW is insufficiently high. After T^ec, the dynamic 
evolution of the shock can be obtained from the energy 
conservation law £ — T 2 M. sw c 2 as 



T(t) 



3£ 



32imm p c 5 t 3 



i s [ £- 3 /8 ; T dcc < t < T ci 

t , T C [ <C t <C Tg , 

i" 3 / 8 , t>T g . 



(9) 

Obviously, the shock deceleration during T oi < t < T g is 
effectively decreased by the energy injection. 

For the afterglow emission of interest, here we only 
concern the shock emission after T cl . Using the above 
dynamic results and following Sari et al. (1998), the 
(isotropically-equivalent) luminosity of the shock syn- 

4 We will test this argument in Section 4 and give some related 
discussion in Section 5. 

5 Here some shock-related quantities are defined as follows: T 
and R = 2T 2 ct are the Lorentz factor and the radius of the shock, 
Msvi = TjirR 3 nrnp is the (isotropically-equivalent) mass of the 
medium swept up by the shock, and n is the number density of the 
circum medium. m v and c are the mass of proton and the speed 
of light. 
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chrotron X-ray afterglows can be given analytically by 

W (P ' 2)/2 



J ~X — / -X.ci x 



T ci < t < T g , 

(a)- <P - 2,/2 (i)-" P - 2V4 , t >T a , 

(10) 

where p is the spectral index of the energy distribution 
of the shock-accelerated electrons. The luminosity at T e i 
reads 

/" h ~in«/« f ^-i (r f-i r \ (3p ~ 2)/4 

(11) 



3. The former LCs have only one break after the initial 
steep decay phase, while the latter LCs have two breaks. 
Using the following smoothed broken and double-broken 
power law functions (w = 3; Liang et al. 2007), 



F x (t)=R 



X.b 



T, 



(i) 



T, 



(i) 



-l/tu 



(13) 



and 



Y /„ c \(2-p)/2,(p-2)/4 



B -2 er S S 

where <7 p = (p — 2) / (j> — 1), i/x is the X-ray frequency, 
and e e and are the equipartition factors of the electron 
internal energy and magnetic energy, respectively. 

Equation shows that, for p being not much higher 
than two, the shock luminosity is mainly sensitive to the 
parameters e e , £, i m( j, and fs- So an upper limit for the 
shock luminosity for t > T<j e c can be given as follows 



v< 2 > 



T, 



(2) 
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-1/v 
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A comparison between Equations ((H) and (TT21 shows that 
the shock emission can exceed the wind emission only 
for e e ^> £/(l — £). By varying the value of £ and fix- 
ing e e = 0.2, we plot the wind- and shock-contributed 
X-ray LCs in Figure 1. Combining the two emission 
components, the model in principle predicts three types 
of LCs as exhibited, i.e., an energized shock-dominated 
type, a magnetar wind-dominated type, and an inter- 
mediate type. In particular, the intermediate type LCs 
could have an interesting but complicated profile. It will 
be a valuable attempt to find some observational cases 
with such a profile, although it would be not easy. For 
a general investigation, we would only concern the other 
two types in this paper. Anyway, it seems soemwhat too 
simple that nearly all of the Swift-XRT afterglows can 
be summarized by a canonical LC. This view obscures 
the possible diversity of the physical origins of the X- 
ray afterglows. In our opinion, we had better treat all 
of the observed afterglows as a collective of several dif- 
ferent classes, which corresponds to different afterglow 
mechanisms. 

4. CONFRONTING THE MODEL WITH 
OBSERVATIONS 

4.1. Fittings to some observational LCs 

As shown by Figure 1, both the energized shock- and 
magnetar wind-dominated afterglow LCs initially have a 
very flat segment, the decay index of which is close to 
be zero if p is not much higher than two. Therefore, we 
select observational samples from the public XRT-team 
LC repository (Evans et al. 2007, 2009) between 2005 
January and 2009 December under four criterions: (i) 
the data are rich enough to clearly exhibit the profile 
of the LC and no flare appears after the initial steep 
decay phase; (ii) a remarkable plateau emission (a < 0.3) 
immediately follows the initial steep decay phase; (iii) the 
decay index of the last segment of the LC is not much 
larger than 2; and (iv) the GRB's redshift is known. As 
a result, 16 representative samples are obtained, which 
naturally fall into two classes as shown in Figures 2 and 



to fit the selected onc-brcak and two-break LCs, respec- 
tively, we can obtain the slopes (a), the break times 
. t and the X-ray fluxes at the first break (Fx,b) of the LCs, 

Co/B,-i-^md,i,47 erg s(12) as listed in Tables 1 and 2. The distributions of these 

fitting parameters are exhibited in Figure 4 (solid his- 
togram), which shows that (i) is statistically a bit 
larger than and (ii) a^ 1 ' could be usually larger 

(2) 

than ~ 1.2 — 1.4, whereas a 2 is inclined to be smaller 
than ~ 1.2 — 1.4, as expected by our model. But when 
we use Equation (fT3"]) to fit only the first two segments 



of the two-break LCs, the difference between and 

(2) 

a 2 becomes ambiguous as shown by the dashed open 
histogram in Figure 4. Anyway, in contrast to the cen- 
tralization of at ~ 10 4 ~ 5 s, the obvious lack at least 
before ~ 10 6 s of the second break in the one-break LCs 
still suggests a clear difference between these two types 
of LCs, even though a longer-term (> 10 6 s) observation 
may be able to detect a second break in some so-called 
one-break afterglows. 

For the above observational results, a plausible expla- 
nation can in principle be provided by our model where 
the one- and two-break afterglows are respectively at- 
tributed to the shock and wind emission, although it 
seems hard to observationally judge (within the present 
small sample) whether there is an intrinsic difference be- 
tween the one- and two-break LCs. On one hand, for 
the temporal decay indices, the values of \ and ol\ \ 3 
are statistically in rough agreement with the model- 
predicted ct\^ and a!i2*3, respectively. On the other 

hand, the similar distributions of and Tj? can also 
be naturally understood by ascribing both of them to 
a same physical origin, i.e., the r-mode instability re- 



quiring T, 



(i) 



T, 



(2) 



61 



(1 + z)T g , Of course, for a more 



comprehensive observational test, an investigation on the 
afterglow spectral information is desired. On the theo- 
retical aspect, it would be useful to derive some closure 
relations between the temporal decay index and the spec- 
tral index in the wind emission model, as did by some 
authors for the shock model (Zhang et al. 2006; Liang et 
al. 2006; Willingale et al. 2007). Unfortunately, such an 
attempt now is restricted by the ignorance of the specific 
radiation mechanisms of the winds. 
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4.2. Implications to the central magnetars 

Following the above comparison between the model 
and the observations, here we simply connect the selected 
one- and two-break afterglows with the shock- and wind- 
dominated emission, respectively, in order to give a sight 
into the properties of the possible existing central mag- 
netars. To be specific, for the two-break afterglows, the 
values of the characteristic timescales in the model can be 
easily determined by T g ~ T^'/(l+z), T c ~ T^ ] /(l + z), 
and T m = (TgTj 1 / 2 (for q = 1). Then the spin periods 
and magnetic field strengths of the magnetars can be 
easily obtained by 



-l/6j.l/6 



3,3 J 



B 14 = 1.4R 



1/2 



(15) 
(16) 



However, for the one-break afterglows, the values of T c 
and thus T m can not be found from the observations. 
So we have to seek help from Equations ([5]), (ITT|) . and 
obs 6 Then we get 



*-'X,ei 



/•obs 6 
'-X.b- 



(JO 



obs \4/(3p-2) 
X,bA6) 



P 



1/2 



(17) 



with model parameters vx,n.5 = 9p, -0.5 — £ e,-i = 
e_B.-2 = £i,50 = /b,-i = 1 and £ ~ 1. Using the above 
equations, we derive the values of the magnetar param- 
eters Pi and B from both the one-break and two-break 
samples, as listed in Tables 1 and 2, respectively. Figure 
5 further shows that, for the model parameters adopted, 
there is no clear separation between the two samples. 
However, keep in mind that the magnetar parameters 
inferred from the one-break afterglows are strongly sen- 
sitive to the uncertain model parameters (especially p, e e 
and f B ). 

So as a conservative treatment, here we only analyze 
the relatively credible magnetar parameters from the 
two-break afterglows. First, the hypothesis of B < B c 
is favored by the inferred magnetic field strengths as 
B ~ 10 14 " 15 G. Secondly, by calculating the magnetic 
dipole luminosity using Equation (|5|) , we can estimate the 
X-ray radiation efficiency of the magnetar winds, which 
reads 



t = f B £fl/L md ,i ~ (0.01 - 0.1)/ B ,_i 



(18) 



as shown in Table 2. Finally, we even find that the mag- 
netic field strengths and the spin periods of the magne- 
tars may satisfy a loose correlation as B oc P^ Z ^ 2 . This 
is qualitatively in agreement with the dynamo origin of 
the magnetic fields (e.g., Xu et al. 2001). 

5. SUMMARY AND DISCUSSION 

Based on two assumptions of that (i) some GRB cen- 
tral objects are millisecond magnetars and (ii) the mag- 
netar winds can continuously produce X-ray emission 



whose luminosity tracks the magnetic dipole luminosity, 
we investigate the temporal behaviors of the GRB X-ray 
afterglows arising from an emitting magnetar wind and 
an energized external shock together. The competition 
between the internal- and external-origin emission com- 
ponents determines the diversity of the observed X-ray 
afterglow LCs. A comparison between the model and 
observations shows that the model-predicted shock- and 
wind-dominated emission is qualitatively consistent with 
the observed one-break and two-break afterglows, respec- 
tively. 

In the conventional shock model, the second break of 
the two-break afterglows is always connected to the jet 
break, which is however seriously challenged by the usu- 
ally observed chromatic breaks (Liang et al. 2008). In 
contrast, the chromatic breaks could be acceptable for 
the internal-origin emission. On one hand, such an argu- 
ment is supported by the lack of the optical counterparts 
of X-ray flares, which are of internal origin. On the other 
hand, in view of the possible small radii where the inter- 
nal dissipations occur, the low-frequency emission of the 
wind is quite likely to be suppressed, for example, by 
some self absorption effects. Therefore, for some GRBs, 
while their X-ray afterglows are contributed by the mag- 
netar winds, the optical emission could be still dominated 
by the external shocks. In this case, chromatic breaks 
would be detected naturally. 

In this paper, we mainly concern the ordinary GRB X- 
ray afterglows that may be associated with a magnetar 
with a relatively low magnetic field (B < B c ). For sim- 
plicity we do not pay much attention to some unusual 
X-ray afterglows such as the afterglows behaving as a 
plateau followed by a steep decay (e.g., GRBs 060607A 
and 070110). As discussed in Section 3, such X-ray af- 
terglows could be internal afterglows produced by the 
magnetars with relatively high magnetic fields (B > B c ). 
In this case, however, if a low value of B is found, the 
related magnetar could be a candidate of strange quark 
stars rather than neutron stars, since only strange stars 
can suppress the r-mode instability effectively (Yu et al. 
2009b). 
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The fitting parameters* of the two-break X-ray afterglow LCs and the corresponding magnetar parameters. 
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Fig. 1. — Illustrative X-ray (yx = 3 X 10 17 Hz) afterglow LCs contributed by the magnetar wind with varying £ as labeled and by the 
energized shock with t e = 0.2 and g Pi _o.5 = _2 = Ei 50 = 1. The magnetar parameters are taken to be Pi = 1.2 ms, B = 4 X 10 14 G, 
and fg = 0.1. 
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Fig. 2.— 



Selected one-break X-ray afterglow LCs, which can be fitted by Equation (1 1 3 1) as shown by the solid lines. 
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Fig. 4. — Distributions of the fitting parameters of th e on e-bre ak ( lower panel) and two-break (upper panel) X-ray afterglow LCs (solid 
histogram), which are respectively fitted b y Eq uations H13H and 114i . The dashed open histogram shows the results of the fitting to the 
first two segments of the LCs by Equation 1131 . 
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Fig. 5. — Distribution in the P; — B panel of the magnetars inferred from both the one-break (open cycle) and two-break (solid cycle) 

— 3/2 

X-ray afterglows. The dashed line represents B14 ~ 7P i _ 3 . 



